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ABSTRACT 

Compton scattering within the accretion column of magnetic cataclysmic variables 
(mCVs) can induce a net polarization in the X-ray emission. We investigate this 
process using Monte Carlo simulations and find that significant polarization can arise 
as a result of the stratified flow structure in the shock-ionized column. We find that 
the degree of linear polarization can reach levels up to ~ 8% for systems with high 
accretion rates and low white-dwarf masses, when viewed at large inclination angles 
with respect to the accretion column axis. These levels are substantially higher than 
previously predicted estimates using an accretion column model with uniform density 
and temperature. We also find that for systems with a relatively low-mass white dwarf 
accreting at a high accretion rate, the polarization properties may be insensitive to 
the magnetic field, since most of the scattering occurs at the base of the accretion 
column where the density structure is determined mainly by bremsstrahlung cooling 
instead of cyclotron cooling. 
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1 INTRODUCTION 

It is anticipated that X-ray polarimetry will provide us with 
a powerful method of probing the physical conditions and 
geometry of high-energy astrophysical systems. Accreting 
X-ray sources are expected to show a significant degree of 
polarization as a result of photon scattering in non-uniform 
distributions of matter in non-spheric al geometries, such as 
accretion disks and columns (see e.g. M eszaros et alJll988t 
lReeslll975h . 

Accreting white dwarfs are strong X-ray sources durin 
active states (see e.g. iKuulkers et al. [ 120061 ; IWu et~ai1l200 
Warner 1995, for reviews). In magnetized systems (the mag- 
netic cataclysmic variables, mCVs), the accretion flow is 
confined by the magnetic field near the white dwarf. The 
supersonic accreting material becomes subsonic close to the 
white-dwarf surface, resulting in a standing shock, which 
ionizes and heats the plasma to temperatures kT « (10 — 40) 
keV, (where k is Boltzmann's constant). The heated mate- 
rial in the post-shock flow cools by emittin g bremsstrahlung 
X-rays and optical/IR cyc lotron emission (|Lamb fe Masters! 
ll979l ; lKing fe LasotdHimSh . Bremsstrahlung radiation emit- 
ted by isotropic thermal electrons is not polarized. How- 
ever, in strongly magnetic systems where cyclotron cool- 
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ing is very efficient, Coulomb collisions might not be effi- 
cient enough to ensure an isotropic Maxwellian distribution 
for the electrons. Bremsstrahlung X-rays from such system s 
would be intrinsically polarized (see e.g. iMcMasteii Il96l[ ) . 
For mCVs with a high accretion rate, the accretion column 
can have Thompson optical depths of up to a few, gi ving 
rise to substantial Comptonization signatures (see e.g. | W ul 
1 19991 ; iKuncic. Wu fc Culled [20051 ; IMcNamara et alj|2008j h 

Our prev i ous studies l|Kuncic et al. 2005; 

IMcNamara et al". 

I 120081 1. which used a nonlinear Monte 
Carlo algorithm for the simulations (|Cullenl l2001al lbh. 
demonstrates the substantial effects of Compton scatter- 
ing on Fe Kq emission lines in the post-shock flows of 
mCVs. Like the photons in the Fe lines, whose profiles 
are broadened and distorted by Compton scattering, 
the photons in the whole X-ray continuum can undergo 
multiple scatterings. Although Compton scattering would 
not readily introduce prominent spectral signatures in the 
0.1 — 10 keV continuum energy band, it can produce a 
net polarization due to the non-isotropic distribution of 
electrons in the accretion column, the lack of symmetry 
in the viewing geometry and perhaps the pre sence of the 
magnetic field in the accretion flow. A study bv lMattl (|2004h 
has shown that the degree of polarization is ~ 4% for a 
cold, homogeneous and static accretion column. The degree 
of X-ray polarization in mCVs should be higher, given that 
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the accretion flow is stratified and has non-zero temperature 
and velocity. It is expected that this polarization will be 
observable. 

In this paper, we investigate the X-ray polarization 
properties of mCVs by means of Monte Carlo simulations. 
We consider a more realistic model for the accretion column, 
which takes into account the full ionization structure of the 
post-shock column. The velocity, temperature and density 
profiles of t he post-shock flow are derived using a model as 
described in IWu. Chanmugam fc Shavivl i| 19941 ) . In Section 
2, we outline the theory of polarized Compton scattering 
and the computational algorithm for the polarization calcu- 
lation. In Sec. 3, we present our findings. A summary and 
conclusion are given in Sec. 4. 



2 POLARIZED COMPTON SCATTERING 
2.1 Basic Physics 

The differential Klein-Nishina cross-section for the scatter- 
ing of polarized electromagnetic wav es off electrons is (see 
lHeitledfl936l ; IJauch fc Rohrlichlll980l . for example), 
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Here, ro is the classical electron radius, -y = (1 — /3 ) _ ' is 
the Lorentz factor of electrons with velocity v = /3c, p is the 
cosine of the angle between the propagations of the incident 
photon and electron, and 



X 



K K 
2 + K 



+2 e-e' + 



- 1 



e ■ pe ■ p 



(2) 



• pe ■ p 



where p is the electron momentum, e is the polarization 
vector of the incident photon and k = —p-k and k' = —p-k' , 
where k is the incident photon momentum. All other primed 
quantities denote values after scattering. Also, 
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where v is the initial frequency of the photon and 6 is the 
scattering angle (angle between the incident and scattered 
photon propagations), and p! is the cosine of the angle be- 
tween the propagations of the scattered photon and the in- 
cident electron. 

The expression for X given by ([2} simplifies consider- 
ably by specializing to the rest frame of the incident electron. 
In this case, p — (m c c 2 , 0,0,0) and a specific gauge can be 
chosen such that e*o = ep ' = 0, then p-e = p • e ' = and 
becomes (see e.g. IJauch fc R ohrlich 1980, for details) 



X ( 



- 1 + 2 (e 



(4) 



where the subscript 'e' denotes a quantity calculated in the 
electron rest frame and 



7/c 



m c c 2 



(1 



(•») 



The differential cross section in the electron rest frame sim- 
plifies to, 
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Figure 1. Schematic geometry of an mCV accretion column with 
total height H and cross-sectional radius r c , viewed from an in- 
clination angle i. 
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We can consider the scattered radiation as composed of 
linearly polarized components perpendicular and parallel to 
the incident plane of polarization. Then the differential cross 
section for polarized scattering in the electron rest frame can 
be written in the more familiar form 
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The polarization vector e" for the polarized fraction P 
of photons is perpendic ular to the scattering plane and is 
defined by (|Angellll969n . 



e' 



(8) 



where Cl' is the propagation directional unit vector of the 
scattered photon. For the remaining 1— P fraction of photons 
e" is randomly distributed in the plane perpendicular to Cl' . 
In the case where the radiation is intially unpolarized, the 
degree of linear pol arization indu ced by Compton scattering 
is given by (see e.g. |Polan|[l96^ . 



P = sin 2 6 C (t? c + % 1 
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2.2 Computational Algorithm 

IPozdnvakov. Sobol fc Sunvaevl (|l983T ) outlined a Monte 
Carlo algorithm to model Compton scattering of photons. 
The algorit hm was i mplem ented into a non-linear code de- 
veloped by ICullenl <|2001ah which was subsequently used 
to investigate Compton scattering of Fe lines in mCVs 
ijKuncic et al.ll2005l : iMcNamara et al.ll200Sr i. Here, we gen- 
eralize the algorithm to include the polarization. If polariza- 
tion is not included, the azimuthal angle <f> of the scattered 
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photons may be sampled uniformly. However, in polariza- 
tion calculations the azimuthal distribution of the scattered 
photons depends on the polarization vector of the incident 
photons and is anisotropic, thus uniform sampling is not 
applicable. In the algorithm for our previous calculations, 
the energy and direction of the scattered photons are calcu- 
lated in the white-dwarf frame. Polarized scattering in the 
white-dwarf frame, however, requires sampling multivariate 
distributions of the scattered photon energy hv' , the cosine 
of the angle between the scattered photon and electron //, 
the azimuthal angle cf> and the angles between the polar- 
ization vectors and photon propagation vectors. To avoid 
this complication, a Lorentz transformation into the inci- 
dent electron rest frame is performed. The Lorentz transfor- 
mation of the photon momentum and energy between the 
electron rest frame and the lab frame is given by 



and 
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u e = vy(l-/3 cos 6) (11) 

respectively. Once the resulting energy and momentum are 
calculated in the electron rest frame, a transformation is 
made into the lab frame using (|10[) with a reversed /3. The 
resulting scattering angle is transformed using 

cos6» c + /3 

cos# = - 7, 7T ■ (12) 

1 + f3 cos 9 C v ' 

In the incident electron rest frame we sample the vari- 
ables fj,' B and <f)a using the inverse-function method and devise 
a rejection tech nique . We follow the same proc edure outlined 
in ISoboll (|l9T9h and IPozdnvakov et all (|l983h to determine 
the joint distribution density of the variables, 

We rewrite the joint density in the form (|Sobollll979l ) 
1104,0,0 =a- 1 Il 1 ( l i' cl( j> c )Y (14) 
where the normalized density Hi(fj,' e , C ) = J- cos 2 <f> e and 
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Our algorithm is as follows: 



(i) In the white-dwarf frame, choose the incident electron 
energy E and momentum p from a Maxwellian distribution. 

(ii) Calculate fi, (3 and 7. 

(iii) Calculate the total cross-section a. 

(iv) Determine whether the scattering event will occur 
using a rejection algorithm (Cullen 2001a, b). If the scatter- 
ing event is rejected, then the photon continues in the same 
direction without scattering. 

(v) If scattering is accepted, then 

(a) Lorentz transform the incident photon momentum 
k, energy hv and \i into the electron rest frame. 

(b) Select two random numbers ui and U2, and cal- 
culate a possible direction of scattering for the density 
fli(/ic, <t> c ) by solving: 



fi c = 2ui — 1 
24> c + sin(20 e ) 



(16) 



(c) Calculate the scattered photon direction Q' e and 
scattering angle 6 C from Cl c ■ Cl' c . 

(d) Determine the energy of the scattered photon hv' e 
from, 



[l + (hv c /m c c 2 )(l - cos e )] 



(17) 



(e) Determine Y from (|15|) . If Y ^ 2, accept the scat- 
tered quantities, else recalculate fi' e and 4> e . 

(f) Once n' e and cf) a are known, determine whether the 
scattered photon will be polarized by calculating the de- 
gree of polarization P and extracting a random number 
uz- If U3 < P, then the photon is unpolarized else the new 
polarization vector e' a is chosen from JS). 

(g) Lorentz transform the scattered photon momentum 
fe'o, energy hv' e and scattering angle B into the white- 
dwarf frame. 

(h) Calculate the total linear polarization P of photons 
emerging in a specified direction as follows: define a view- 
ing plane and calculate the Sto kes parameters Q and U 
by projecting e onto the plane (|Matt et al.l 1199a ). Then 
calculate P = ■ S /Q 2 + U 2 /I, where I is the number of 
photons scattered in the chosen direction. 

In this paper, we only consider unpolarized incident 
photons and the initial polarization vector e is chosen ran- 
domly in the plane perpendicular to the photon propagation 
direction Cl, i.e. with the condition e • = 0. The energy of 
the incident photons is sampled from a bremsstrahlung spec- 
tral function N(v )dv using rejection methods as described in 
iPress et al](|l992l ). Since the distribution function N{v)dv is 
too complex to sample directly, we consider an approximate 
spectral function F(v)dv which lies above N{v)du every- 
where for a predefined energy range /w max ^ hv ^ /w m j n . 
We use the comparison function from lFromerth et alj (|200ll ) 
for the thermal bremsstrahlung photon distribution function 
with cutoff temperature T, 



12 for hu mia /kT < x < 1, 
for x > 1. 
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where x = hv/kT. The comparison function is sampled to 
determine the energy of a photon hv. If hv lies in the prede- 
fined energy range and uf(y) ^ N(v), where u is a standard 
deviate, then the sampled energy is accepted. If the condi- 
tions are not satisfied then we sample a new energy from the 
comparison function. 



3 RESULTS AND DISCUSSION 

Our model scattering region is a cylindrical column, which 
is divided into a heated post-shock flow and a cool pre- 
shock flow. Since the upstream infall is supersonic and close 
to free fall, we assume that the velocity and density are 
constant in the pre-shock region. The post-shock flow is 
stratified in density, velocity and temperature, and these 
profiles ar e calculated usin g the hydrodynamic model de- 
scribed in IWu et all l| 19941 ). The relative efficiency of cy- 
clotron and bremsstrahlung cooling is determined by the 
magnetic field, through an efficiency parameter e s evaluated 
at the shock. The ionization structure of the post-shock re- 
gion is dependent on the white-dwarf mass Mwd, and ra- 
dius Rwd, the specific mass accretion rate rh and the ratio 
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Photons emitted Photons initially propagating across the Photons which subsequently escape 

isotropically near column have a higher probability of escaping, towards i = 90 do so as a result of large 

column base. since the optical depth is lower. angle scatterings, as the photon distribution 

becomes increasingly anisotropic. 



Figure 4. Propagation of an isotropic distribution of photons injected at the base of the mCV accretion column. As the photons 
propagate through the length of the column their distribution becomes increasingly anisotropic. 
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Figure 2. Polarization degree P plotted as a function of the in- 
clination angle i for an accretion column with a white-dwarf mass 
■^WD = l.OAfg, specific accretion rates of rh = l.Ogcm - 2 s _1 
(dotted line) and m = 10gcm~ 2 s _1 (solid line) and for different 
column hcight-to-radius ratios H/r c , as shown. The optical depth 
across the column at the shock is r = 0.3 and 0.04 for the high 
and low rh cases, respectively. 



Figure 3. Polarization degree P plotted as a function of the in- 
clination angle i for an accretion column with a white-dwarf mass 
Mwb = 0.5 Mq, specific accretion rates of rh = 1.0gcm _2 s _1 
(dotted line) and m = 10gcm _2 s _1 (solid line) and for different 
column hcight-to-radius ratios H/r c , as shown. The optical depth 
across the column at the shock is r = 1.0 and 0.1 for the high 
and low rh cases, respectively. 
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Figure 5. Same as Fig. 02(b), except for a uniform density ac- 
cretion column 



Figure 8. Linear polarization P plotted as a function of phase 
angle <t> for the mCV GK Per with a white-dwarf mass Mwd = 
0.63 Mq and specific accretion rate rh = 10 gem - 2 s — . 
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Figure 6. Polarization degree P plotted as a function of inclina- 
tion angle i for an mCV with a white-dwarf mass A/wd = 0.5 Mq 
and specific accretion rate rh = 10 gem -2 s — 1 with a ratio of cy- 
clotron to bremsstrahlung cooling e s = 10 (solid curve) and e B = 
(dotted curve). 
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Figure 7. The average number of scatterings per photon plot- 
ted as a function of inclination angle i for a white-dwarf mass 
Mwd = 0.5 Mq with H/r c = 5. The solid and dotted line shows 
the average scattering per photon for rh = 10 and lgcm~ 2 s _1 , 
respectively (c.f. Fig. Qjjj). The dashed line shows the results for 
rh = 10gcm _2 s — 1 when e s = 10 (c.f. Fig. [BJ. The dot-dash line 
shows the average number of scattering per photon for a uniform 
column with rh = 10 g cm -2 s — 1 (c.f. Fig. [5]l . 



of t he efficiencies of cyc l otron to bremsstrahlung cooling e s 
(see iKuncic et ail 120051 ; iMcNamara et all I2008I . for futher 
details). In this study, we consider cases where e s = and 
10. We also assume that the emitted bremsstrahlung pho- 
tons are intially unpolarized and are emitt ed from t he base 
of the column where the emissivity peaks l|Wull 19941 ). 

We consider an accretion column for two different white- 
dwarf mass- radius values: Mwd = 0.5 Mq, Rwd = 9.2 x 10 s 
cm a nd Mwd = 1.0 Mq, Rwd = 5.5 x 10 8 cm (|Nauenberd 
Il972l) . For each mass, we consider two different specific mass 
accretion rates: rh = lgcm _2 s _1 and rh = 10gcm _2 s _1 , 
roughly corresponding to the low and high ends of typical 
accretion rates of mCVs. The viewing inclination angle i is 
measured from the cylinder axis (see Figure [1J. A polariza- 
tion plane is defined for a small range of inclination angle 
i such that Ai ~ 5°. In the simulations, each photon is fol- 
lowed until it escapes the column. If a photon strikes the 
polarization plane, it is binned and contributes to the over- 
all polarization of the emerging beam for that particular i. 

The degree of polarization is very sensitive to the ge- 
ometry and density structure in the accretion column. For 
a lower Mwd, corresponding to a larger Rwd, the optical 
depth across the column is higher for a fixed rh and thus, 
photons undergo more scatterings on average. In the case 
of multiple scatterings, the resulting polarization is largely 
determined by the conditions at the last scattering surface. 
We consider a range of different column geometries by vary- 
ing the column height, H. The column radius r c is fixed to 
0.1-Rwd. 

Figures [5] and show the linear polarization P plotted 
as a function of i for white-dwarf masses Mwd = 1.0 Mq 
and Mwd = 0.5 Mq, respectively with e s = 0. The solid 
curves are for rh = 10 gem -2 s _1 and the dotted curves are 
for rh — lgcm _2 s _1 . We consider three different column 
geometries: H/r c = 3, 5 and 10. In each case, 10 9 photons 
were injected and emitted isotropically at the base of the 
accretion column where the density and hence, emissivity 
peaks. We define P as positive when the polarization is per- 
pendicular to the projection of the cylinder axis onto the 
projection plane, as is the convention for axisymmetric ge- 
ometries. 

In all cases, P increases with rh. This is because the 
electron number density and hence, scattering probability 
increases with rh. Indeed, in all the low rh cases, the average 
number of scatterings is 1, whereas photons scatter on av- 
erage more than once in all the high rh cases (see Fig. 0. The 
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scattering probability is also slightly higher for larger H/r c , 
producing higher P for these column geometries. Figs. [2] 
and [3] also reveal that P increases with i. This trend arises 
because photons in the initial seed distribution at the col- 
umn base escape quickly with none or very few scatterings 
across the radius of the column (z « 90° direction), leav- 
ing behind a photon distribution that becomes increasingly 
more anisotropic (see Fig. |4)| . Photons that now escape in 
directions i m 90° must do so through large-angle scatter- 
ings which result in a high net P. Conversely, photons do 
not need to undergo large-angle scatterings in order to es- 
cape in directions i ~ 0, hence P is lowest for these viewing 
angles. 

The steep rise in P towards i m 90°, reaching values 
up to 8% in the high m cases, differs considerably from the 
results of lMatd (|2004h , which show a gradual rise in P with 
i towards a maximum of 4% at i ~ 90°. This difference 
can be attributed to the assumption of a uniform column 
in Matt's model. To demostrate this difference, in Fig. [5]we 
plot P vs. i for a uniform column with the same parame- 
ters as those used for Fig. [3]d. In this case, the photons are 
emitted uniformly throughout the column. The resulting P 
at i f» 90° is substantially lower than that predicted for the 
nonuniform column (where the photons are emitted at the 
base of the column) . This is because the photon distribution 
remains quasi-isotropic throughout the uniform column and 
thus, fewer photons undergo large angle scatterings to escape 
in the direction i w 90°. Fig. [7] shows the average number of 
scatterings each photon undergoes before escaping the col- 
umn at a particular i. In the case of a uniform column the 
average number of scatterings is lower than the non-uniform 
case for the the same m. 

We also investigate the effect of cyclotron cooling dom- 
inated accretion flows on polarization. Fig. [5] compares the 
polarization degrees for an mCV with Mwd = 0.5 Mq, 
rh = 10gcm -2 s -1 and H/r c = 5 when cyclotron cooling 
at the shock is negligible (e s = 0, dotted curve) and when 
it dominates (e s = 10, solid curve). Interestingly, there is 
not much difference between these two cases. This can be 
understood as follows. The presence of cyclotron cooling 
makes the accretion column more compact and the den- 
sity enhancement increases the average number of photon 
scatterings (Fig. 0. However, for a sufficiently high accre- 
tion rate, the column will be optically thick regardless of 
whether bremsstrahlung cooling or cyclotron cooling dom- 
inates at the shock. For multiple scatterings, the resulting 
polarization degree is determined largely by the final scat- 
tering. Generally in cyclotron dominated flows, the last scat- 
tering surface will be in the dense bremsstrahlung cooling 
zone close to the base of the column, where the physical 
conditions are very similar to the bremsstrahlung cooling 
dominated case. In both cases, the scattering electrons are 
relatively cool with low velocities, and as a result aberra- 
tion is negligible and Compton recoil will dominate. With 
three or four scatterings beforehand, incident photon direc- 
tions for a subsequent final scattering will be approximately 
randomized and hence unpolarized. Thus, the average po- 
larization P of photons emerging from the last scattering 
surface depends only on the scattering angle (and hence the 
viewing inclination) and is insensitive to cyclotron cooling. 

The situation may be different for lower accretion rates. 
If the post-shock region becomes semi-transparent to scat- 



tering, photons can emerge from the entire structure of the 
post-shock flow. In that case, cyclotron cooling, and by im- 
plication the white-dwarf magnetic field which modifies the 
flow, can affect the polarization. We may thus conclude that 
polarization is a robust diagnostic of geometry and orienta- 
tion of systems if they have sufficiently high accretion rates. 

Finally, we show an application to an mCV with high 
accretion rate and low magnetic field. We construct a 
model accretion column using mCV parameters similar to 
that of the intermediate polar GK Per (Fig. [51): Mwd = 
0.63 M(7) and m = 10gc m~ 2 s~ 1 (|Morales Rueda et al.l 
120021 ; IVrielmann et all 120051 ). The predicted fractional po- 
larization increases towards $ = 0.33 where P peaks at 
« 0.075. Whether the predicted P and its phase-dependent 
variation could be observed would depend upon, among 
other things, the extent of dilution by background emission. 
Note that the predicted P for GK Per type mCVs would 
be higher as the accretion flow resembles a curtain rather 
than a cylinder. Also, X-rays reflected by dense material 
further upstream in the accretion disk and the disk mag- 
netosphere coupling region will contribute to the observed 
polarization. Note that for very strong field systems with 
low accretion rates (e.g. polars in an intermediate accretion 
state), an anisotropic temperature distribution due to ther- 
mal decoupli ng between charged particles at the shock (see 
ISaxtonll2005h may complicate the situation by imprinting an 
intri nsic polarizati on on the seed bremsstrahlung radiation, 
(see IWarnerlll995l . for a review of polars and intermediate 
polars). Also, an additional polarized Compton component 
resulting from reflection off the white dwarf surface can con- 
tribute to the overall spectrum above a few keV and could 
intro duce an energy dependence on the polarization (|Mat j 
l2004h . 



4 SUMMARY AND CONCLUSION 

We have investigated the properties of Compton polarized 
X-rays in the accretion column of mCVs using Monte Carlo 
simulations. The calculations take into account a post-shock 
region stratified in temperature, density and velocity. The 
degree of linear polarization for scattered X-rays was calcu- 
lated for a range of different column geometries, white dwarf 
masses and accretion rates. 

We have found that the resulting polarization is sensi- 
tive to the density structure in addition to the viewing ge- 
ometry. The non-uniform density structure in the post-shock 
column has a significant effect on the photon distribution 
and average number of scatterings. We have demonstrated 
that enhanced emissivity near the base of the column re- 
sults in photon distributions that become increasingly more 
anisotropic throughout the column length so that X-rays es- 
caping at large angles with respect to the column axis can do 
so only through increasingly larger angle scatterings. This 
produces more strongly polarized X-rays emerging from di- 
rections perpendicular to the column axis than directions 
parallel to it. 

We find that the degree of polarization is two times 
larger than the levels ~ 4% predicted in previous stud- 
ies, for a cold, static accretion column with uniform den- 
sity. We demonstrate that for usual mCV parameters (e.g. 
Mwd = 0.5Mq and m = 10gcm _2 s _1 ), the degree of 
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polarization can reach up to ~ 8% for accretion columns This paper has been typeset from a TjjX/ P/TffX file prepared 
with high viewing inclination angles. The polarization is not by the author, 
significantly affected by introducing cyclotron cooling, pro- 
vided that the accretion rate of the system is sufficiently 
high. 
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